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Abstract 

The response of an array of plastic phoswich detectors to ions of 1 < Z < 18 
has been measured from E/A=12 to 72 MeV. The detector response has been 
parameterized by a three parameter fit which includes both quenching and 
high energy delta-ray effects. The fits have a mean variation of < 4% with 
respect to the data. 



I. INTRODUCTION 



In recent years a number of large arrays using plastic and/or Csl scintillators have been 
built for heavy ion experiments |l -16]. While these arrays offer a relatively inexpensive way 
of covering large solid angles, the nonlinear response of various scintillators as a function of 
the incident particle's energy, charge, and mass poses a challenge in their calibration. To 
understand the nonlinear response of various scintillators, calibration data covering a wide 
range of projectile energy, charge and mass are required. These requirements, combined 
with the large number of elements in an array, make it virtually impossible to do a detailed 
calibration after each experiment. 

Thus a procedure by which the calibration of large scintillator arrays can be achieved 
during or after an experiment with only a few calibration points would be valuable. The 
essential ingredient in such a procedure is an understanding of the scintillator response to a 
wide variety and energy of incident particles. Considerable work has been done in the past 
to understand the complex scintillator response. Birks proposed [17]] a simple expression for 
the light output which took into account the quenching of light output for cases of very large 
energy loss. Other authors proposed expressions which also took into account high energy 



electrons, or "delta rays", which escape from the primary ionization column ]p!8|-p0|, and 
recently a new model of light production by secondary electrons has been suggested |2"T| , |2"2" 



While some authors have applied these models to their data p3| , ^4| , many others have used 
empirical expressions which may not be generally applicable [^,|6iP|, [lQ|,p5| -|30[ . 

In the present work we extend our previously published results for Csl response |U| by 
incorporating delta rays into the Birks' formalism and applying the resulting expression to 
calibration data taken with plastic phoswich detectors from the Chalk River/Laval forward 
array ||. Details of the calibration data are given in the next section. In Section |T| the 
parameterization of the light output is derived. Results with fits to the parameterization of 
Section |T| to the calibration data are discussed in Section [TV]. A step by step description of 
the procedure used to apply the calibration to data is described in Section [V] Conclusions 
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from the present work are given in Section 



VI. 



II. MEASUREMENTS 

The CRL/Laval forward array phoswich detectors || consist of 0.7 mm thick fast plastic 
AE (BC408 from Bicron Corporation) elements heat-pressed to 76 mm thick slow plastic 
E (BC444) elements. The photomultiplier signal is split, and integrated with two charge- 
to-digital converters (QDCs). The first QDC, with a gate width of 57 ns, integrates the 
early portion of the signal which corresponds primarily to energy lost in the fast plastic AE 
counter. The second QDC, with a gate width of 183 ns is delayed by 210 ns with respect to 
the first gate. This QDC integrates a later portion of the signal which corresponds primarily 
to the energy deposited in the slow plastic E counter. 

An example of a AE — E particle identification plot is shown in Fig. [I]. Ions which stop 
in the fast plastic AE lie along the steeply rising line, or "backbone", near the y-axis. The 
deviation of the backbone from the vertical is due to the tail of the fast plastic signal which 
extends into the delayed slow gate. Neutral particles, almost all of which interact only in 
the slow plastic E, are visible in the inset in Fig. [I] as the straight line below the hydrogen 
isotopes. The nonzero slope of the neutral line is due to the leading edge of the slow signal 
which lies within the fast gate. 

In order to measure the response of phoswich detectors from the CRL/Laval forward 
array, four elements of the array were mounted on a movable trolley which allowed them to 
be placed directly in the beam. The detectors were exposed to a series of secondary beams 
produced from an 18 beam at E/A=3Q.5 MeV incident on two 12 C production targets 
at the exit of the cyclotron at the TASCC facility at Chalk River. Six beamline rigidity 
settings, listed in Table [I], were used to select different energies and isotopes. The production 
target-beamline rigidity combinations produced up to 94 calibration points per detector. A 
total of 41 different isotopes with 1 < Z < 10 and 12.0 < Ej A < 72.6 MeV were produced. 
Representative AE — E plots are shown in Fig. |2] for two beam rigidities. Only calibration 
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points with more than ten counts were used in the subsequent fitting procedure. 

In order to convert the centroids of the calibration points into values proportional to the 
light emitted in each element of the phoswich, the QDC pedestals must be subtracted, and 
a correction must be applied for the cross talk of the AE and E signals in the two QDC 
gates. The correction for the latter signal can be determined from the slopes and intercepts 
of the backbone and neutral lines. We take 

s = L + m b dL + b b (1) 

and 

f = dL + m n L + b n , (2) 

where s and / are the measured E (slow) and AE (fast) signals, L and dL are the "true" 
E and AE signals, m n (m b ) is the slope (inverse slope) of the neutral line (backbone), and 
b b and b n are the QDC pedestals. Solving for dL the measured signals can now be related 
to L by: 

L _ s - m b f + m b b n - b b 
1 - m b m n 

In this paper we describe only the calibration of the slow plastic signal. The total energy 
may then be obtained by determining the energy loss in the AE element from a table lookup 
based on the ion's element number and the energy deposited in the slow plastic. 

III. LIGHT OUTPUT PARAMETERIZATION 

In order to account for the observed nonlinearity of scintillator light output Birks pro- 
posed a simple expression for the differential light output ^ in terms of a quenching coeffi- 



cient kB times the energy loss j^: 



dL S dE 

dx ~ 1 + ( ' 
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where S is the scintillation constant [[17]]. Under this formulation, an increasing energy loss 
leads to greater quenching as the primary ionization column becomes saturated. The next 
step is to consider the effect of high energy electrons, or "delta rays" , which have sufficient 
energy to escape from the primary ionization column. These delta rays, whose scintillation 
efficiency is nearly 100% [|I^], will begin to dominate the light output for very heavy ions 
where the light output within the primary ionization column is almost totally quenched. If 
F s is the fraction of the energy carried by delta rays then the differential light output may 
be expressed as 

dL _ S(l - F s ) dE dE 

dx 1 + (1 - F s )kB^ dx 3 dx' U 
The fraction F s may be expressed as [15 



= l ln(2m e c 2 7 2 /? 2 /T )-/3 2 
s 2 ln(2m e c 2 7 2 /? 2 //) -/3 2 ' 1 ' 

where (3 and 7 are calculated from the ion velocity, m e is the electron rest mass, To is the 
minimum electron energy needed to escape from the primary ionization column, and / is the 
ionization potential of the scintillator, I ^0.048 keV for plastic scintillator. At intermediate 
energies 7 2 ~1, (3 2 <C 1, and Eq. (^) simplifies to: 
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In hE/A\ 



where a = and m is the nucleon rest mass. Since F s is negative for E/A < Eq. 



m 

t. 



is integrated up to E/A = with F s = 0. The light output may now be expressed as: 

^ f ^ okB f* (l + i?) 2 f 



L = SE- SkB / a — — -TpdE - S — / ' d \,~ dE (8) 

Jo 1 + 2 J^2 + kB(l + R)!§ KJ 

where R = t T^Tj^n ' ^' iere are three parameters in this expression: the scintillation constant 
S, the quenching factor kB, and the electron kinetic energy cutoff T . The quantity S is 
measured in combination with an overall gain factor from the readout device, but the latter 
two parameters are, in principle, properties of the scintillator material and should be the 
same for each detector. 



IV. RESULTS 
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O Calibration 



The free parameters in Eq. (Bp were fitted to the data described in Section [TJ for the 
energy deposited in the slow plastic. The energy lost in the fast plastic was calculated 
using the energy-loss code STOPX |33| and subtracted from the incident energy to give the 
energy deposited in the slow plastic. In order to obtain satisfactory fits it was necessary 
to separately fit the light, 1 < Z < 3, and heavy, 4 < Z < 10, ions. The values for the 
electron kinetic energy cutoff parameter T , from the four detectors, were found to be in 
good agreement for the light ion fit, see Table 0. The fit was then repeated with To fixed 
at the average value of 2.85 keV. The results of this second fit showed a better agreement 
in the four values of kB without significantly affecting the quality of the fit. A final pass 
was made with kB at the average value of 8.25. Again, there was no significant effect on 
the overall quality of the fit. The mean difference between the fit and data is < 4%. The 
light ion data and final fit results are shown in Fig. [| for one detector. 

In fitting the heavy ion data, 4 < Z < 10, to Eq. (f|) it was necessary to add a quadratic 
term to Eq. (H) to account for a slight curvature of the backbone for high AE. Equation |3| 
becomes 

s - m b f + m b b n - b b 2 

L = : qf ■ 9 

1 — m b m n 

The additional parameter q may be included in the fit, or determined independently by 
making the backbone pass through a punch-in point of a high Z line such as the one indicated 
in Fig. |l|. Once again, the values for To and kB were found to be in good agreement from 
detector to detector, see Table [[II]. The mean difference between the fit and data is < 2%. 
The results of the fits for the heavy-ion data are shown Figure ^ for one detector. 

After fitting both the light- and heavy-ion data it was found that the ratio of the light 
ion gain factor Sli to the heavy ion gain factor Shi was Ril.30 for all four detectors; see 



Table |TV|- The ability to obtain a single detector-independent set of parameters kB and Tq 



combined with a fixed ratio for the overall gain parameters, permits the calibration of other 
data sets with a single calibration point. 



B. Gate Tests 

The sensitivity of the gains to changes in the slow gate width and delay were checked by 
separately varying the gate width and delay from their normal values. Figure g|(a) shows the 
percentage shift when the gate width was changed from the usual value of 183 ns. Reducing 
the gate width to 150 ns resulted in an average shift of ~ —13.5%, with no apparent Z 
dependence and less than a 1% scatter in the individual shifts. Lengthening the gate width 
to 210 ns produced an average shift of ~ +11.5% with a small Z dependence. The effect 
of varying the slow gate delay from the normal value of 210 ns is shown in Fig. 0(b). For 
changes of ±30 ns, very small (< ±1%) Z dependence was observed. Reducing the slow 
gate delay to 150 ns produced the largest observed Z dependence. Even in this case the 
observed Z dependence is still small, < ±2.5%, with the largest shifts observed for Li and 
Be isotopes. 

Data taken with phoswich detectors can be affected by differences in timing between 
detectors. Since timing differences result in a small smear of slow gate delays, the data in 
Fig. |5|(b) may be used to estimate the effect on the measured energy of poor timing between 
phoswich detectors. In the past we have experienced timing differences of up to 10 ns for a 
few of the phoswich detectors which has now been corrected. From Fig. 0(b) the effect of a 
10 ns time shift can be estimated to be less than 5%. 



C. 37 C1 Calibration 

The extent to which the values of kB, T , and Slj/Shi, which were obtained in Sec- 
tion [IV A|, are applicable to other data sets was checked by calibrating a second data set 



with the parameter values obtained from the first data set. In the second data set, a single 
phoswich from the first ring of the array was placed directly into a secondary beam ranging 



from Hydrogen to Argon ions produced by a 37 CI beam incident on a 50 mg/cm 2 C produc- 
tion target with a beamline Bp setting of 1.584 Tm. With kB and T fixed at 7.18 and 1.13 
keV respectively, the data points for 4 < Z < 18 were fit with just two free parameters, 
Shi and q. The fit results differed from the measured values by an average of 1.2%. The 
light ion data were then fit with kB=8.25 and T =2.85 keV and with the light-ion gain 
factor, Sli, fixed at 1.30Shi- The average difference between the predicted and measured 
centroids for 1 < Z < 3 was 5.6%. While somewhat better fits were obtained when all the 
fit parameters were allowed to vary, the calibrations obtained with the 18 fit parameters 
are quite satisfactory for plastic phoswich detectors. 

V. APPLICATION TO DATA 

The parameterization described above in Eq. [8] is for the energy deposited in the slow 
plastic E detector. In order to obtain the total incident energy from this parameterization it 
is necessary to calculate the energy deposited in the fast plastic AE detector. Furthermore, 
Eq. [8] can only be numerically inverted into a relation for energy as a function of light. The 
following steps are followed to produce a table of the total incident energy as a function of 
the light output of the slow plastic for a given ion: 

1. The energy loss SE in the fast plastic AE detector is calculated for the incident 
fragment energy E inc . 

2. Equation || is used to calculate the light output for the energy deposited in the slow 
plastic E s i ow = E inc - 5E. 

3. Steps 1 and 2 are repeated for a range of fragment energies. The results are used to 
build a table of light output as a function of incident energy. 

4. The table in step 3 is inverted to give a table of the incident energy as a function of 
the light output of the slow plastic. 

5. Steps 1 to 4 are repeated for each ion to which the calibration will be applied. 
Since isotopic resolution in the phoswichs is only acheived in the calibration data through 
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Bp separation, normal data sets use calibrations assuming A = 2Z for Z > 2 and A = 1 for 
Z — 1. This assumption introduces an error in the fragment energy when the fragment has 
a different mass number than the assumed value. For example the calculated energy will 
be «10% less than the true energy for 8 Li and «5% less for 18 0. For comparison the slow 
plastic energy resolution is «5%. 

VI. CONCLUSIONS 

The response of plastic phoswich detectors from the CRL/Laval forward array has been 
measured for a wide range of incident particles and energies. The data are well reproduced 
by a fit based on the light output relation of Birks with an additional term to include the 
effects of high energy delta rays which escape from the primary ionization column. The light 
output is characterized with just three parameters: the gain factor S, the quenching factor 
kB, and an electron cut-off energy T . The quenching factor and electron cut-off energy 
have been found to be constant from detector to detector. In order to achieve reasonable 
fits over a wide range of ions, 1 < Z < 18, it was necessary to fit the light (Z < 3) and 
heavy (Z > 4) ions separately. The inability to fit all ions with the same set of parameters 
may be due to pulse shape differences between the ions combined with the integration of 
only part of the photomultiplier signal. The mean difference between the fits and the data 
is ~4% for light ions and ss2% for heavy ions. The ratio of the gain factors for the separate 
light- and heavy-ion fits has been found to be a constant. It is thus possible to obtain 
future calibrations with a single normalization point. Tests demonstrated that changes in 
the QDC gates have little effect on the calibrations other than in changing the overall gain 
factor. The calibration procedure outlined in this paper has also been successfully applied 
to Csl detectors over a wide range of incident particles from Z=2 to 32. 
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FIGURES 

FIG. 1. AE (fast signal)-^ (slow signal) plot for 70 Ge+ na *Ti at E/A = 35 MeV at lah = 13.4°. 
The inset shows an enlargement of the lower left corner. 

FIG. 2. AE (fast signal)—^ (slow signal) plots for secondary beams for (a) Bp = 1.50Tm and 
(b) Bp = 2.Y7Tm. Selected isotopes are indicated. 

FIG. 3. Final results of fits of Eq. (8) to the light ion data for detector #34. The data are 
shown as points, the results of the fits are indicated by lines. 

FIG. 4. Selected fits of Eq. (8) to the heavy ion data for detector #34. The data are shown as 
points, the fit results are shown as lines. 

FIG. 5. Sensitivity of the detector gains to changes in the slow gate width and delay. The data 
are from one detector for runs with Bp=1.5 Tm. Multiple isotopes are included for some fragment 
charges, (a) Percentage shift in the detector gains for slow gate widths of 150 and 215 ns relative 
to the normal slow gate width of 183 ns. (b) Percentage shift in the detector gains for slow gate 
delays of 150, 178, and 240 ns relative to the normal slow gate delay of 210 ns. 
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TABLES 

TABLE I. Beamline rigidity, Bp, 12 C production target thickness, charge range of measured 
fragments, and degraded 18 energy. For Bp=2.17 Tm the energy for fragments with Z/A=0.5 is 
given. 



Up 


Production Target 


Z range 


yj r^nergy 




( TT\ CT / CTYi ^ ^ 




I A/TpV / nnplpn'r) i 


1.25 


307 


1-9 


14.8 


1.33 


307 


8 


16.6 


1.50 


80 


1-10 


21.2 


1.75 


80 


2-8 


28.8 


1.86 


80 


8 


32.3 


2.17 


80 


1-5 


25.0 (Z/A=0.5) 
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TABLE II. Light ion (1 < Z < 3) fit parameters To and kB, and average percentage difference 
between the fit and the slow-plastic data for each detector. The average percentage difference 
includes only data points with 10 or more counts which are at least 40 channels beyond the 
backbone. The electron cut-off energy To was fixed after the first pass, and the quenching coefficient 
kB was fixed after the second pass. 







Pass 1 




TV...,, O 

_Pass 2 






rass 6 




Detector 


T 


kB 


Difference 


T 


kB 


Difference 


T 


kB 


Difference 




(keV) 




(%) 


(keV) 
(fixed) 




(%) 


(keV) 
(fixed) 


(fixed) 


(%) 


33 


2.84 


8.13 


3.9 


2.85 


8.11 


3.9 


2.85 


8.25 


3.9 


34 


3.24 


7.36 


3.0 


2.85 


7.89 
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2.85 


8.25 


2.9 


35 


2.71 


8.27 
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8.04 
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2.85 
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3.3 


48 


2.60 


9.46 


2.2 


2.85 


8.94 


2.4 


2.85 


8.25 


2.4 


Average 


2.85 


8.31 




8.25 
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TABLE III. Heavy ion (4 < Z < 10) fit parameters, To and kB, and average percentage 
difference between the fit and the slow-plastic data for each detector. The average percentage 
difference includes only data points with 10 or more counts which are at least 40 channels beyond 
the backbone. The electron cut-off energy Tq was fixed after the first pass, and the quenching 
coefficient kB was fixed after the second pass. 
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(keV) 
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(keV) 
(fixed) 




(%) 


(keV) 
(fixed) 


(fixed) 


(%) 
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1.02 


7.18 
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7.18 
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6.44 


1.6 


1.13 


6.43 


1.6 


1.13 


7.18 


1.7 


Average 


1.13 


7.18 




7.18 





TABLE IV. Gain factors from light-ion, Sli, and heavy- ion, Shi, fits to Eq. (8), and the ratio 
Sli/Shi- 



Detector Sli S H i S L i/S H i 

33 4.54 3.52 1.29 

34 4.26 3.33 1.28 

35 5.11 3.94 1.30 
48 4.69 3.58 1.31 

Average 1.30 
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